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Abstract 



We present the results of a search for pair production of scalar top quarks (t\) 
in an i?-parity violating super symmetry scenario in 106 pb _1 of pp collisions 
at y/s =1.8 TeV collected by the Collider Detector at Fermilab. In this mode 
each t\ decays into a r lepton and a b quark. We search for events with two 
r's, one decaying leptonically (e or fj,) and one decaying hadronically, and 
two jets. No candidate events pass our final selection criteria. We set a 95% 
confidence level lower limit on the t\ mass at 122 GeV/c 2 for Br(ii — ► rb) = 1. 
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Many supersymmetry (SUSY) models [1] predict that the first two generations of super- 
symmetric partners of the quarks and the leptons (squarks and sleptons) are approximately 
mass degenerate. However, the mass of the lightest top squark (ti or 'stop') can be relatively 
light due to a large mixing between the interaction eigenstates, II and Ir. This mixing de- 
pends in part on the top Yukawa coupling which is largely due to the heavy top quark mass, 
and it is possible that t\ is lighter than the top quark [2]. 

imparity (R p ) is a multiplicative quantum number defined as R p = (_i) 3B + L + 2S, 5 where 
S, B and L are the spin, baryon and lepton numbers of a particle, respectively [3]. R p 
distinguishes SM particles (R p = +1) from SUSY particles (R p = —1). Conservation of R p 
requires SUSY particles to be produced in pairs and to decay, through a cascade, to SM 
particles plus the stable lightest supersymmetric particle. The Rp conservation, which is not 
required by SUSY, is often built into the theory by hand and is justified phenomenologically 
by limits on the proton lifetime, the absence of flavor-changing neutral currents, etc. Viable 
R p violating [ftp) models can be built by adding explicit ftp terms with trilinear couplings 
(Xijk, Kjki Kjk) an d spontaneous ftp terms with bilinear couplings (e,) to the SUSY La- 
grangian [4,5], where i, j and k are the generation indices. These couplings allow B or L 
violating interactions and, if \' zzk or e% is non-zero, a t\ may decay directly to SM final states 
which are experimentally observable. 

At the Fermilab Tevatron, in pp collisions, stop pairs might be produced strongly via 
.Rp-conserving processes through gg fusion and qq annihilation. In ftp scenarios each stop 
can decay into a tau (r) lepton and a bottom (6) quark with a branching ratio, Br, which 
depends on the coupling constants of the particular model. A good final state search topology 
identifies either an electron or a muon (£ = e or /i) from the r — > tviv r decay, as well as a 
hadronically decaying tau {jh) lepton, and two or more jets. 

We present the results of a search for titi — > £r h jj events, in the framework of Ijtp- 
MSSM, using 106 pb _1 of pp collisions at y/s = 1.8 TeV collected by the Collider Detector 
at Fermilab (CDF) during the 1992—95 run of the Tevatron (Run I). CDF is a general 
purpose detector and has been described in detail elsewhere [6,7]. We briefly describe the 
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subsystems of the CDF detector relevant to this analysis. The location of the pp collision 
event vertex (z vtx ) is measured along the beam direction [8] with a time projection chamber. 
The transverse momentum (p T ) of charged particles is measured in the region \rj\ < 1.0 
by a central tracking chamber (CTC) which is immersed in a uniform 1.4 T solenoidal 
magnetic field [8]. Electromagnetic (EM) and hadronic (HAD) calorimeters, segmented in 
a projective tower geometry surrounding the solenoid and covering the region |^| < 4.2, 
are used for identification of electrons, taus, and jets and the measurement of the missing 
transverse energy (fir). The central strip chamber (CES) is embedded in the central EM 
calorimeter at a depth of approximately shower maximum, and is used for further electron 
identification as well as tt° — > 77 identification from 77 decays. A muon subsystem is located 
outside the hadron calorimeter and has trigger coverage for the region \r)\ < 0.6. 

The analysis begins with a sample of events which pass a three-level trigger system [6] 
which requires a single isolated lepton (e or p) with p T > 8 GeV/c (\r]\ < 1.0 if it is 
an electron and \r]\ < 0.6 if it is a muon) [9]. Offline, the lepton is required to have 
Pt > 10 GeV/c, come from the event vertex, and pass more restrictive identification and 
isolation requirements [7,10]. An event is removed as a Z boson candidate if it contains a 
second, loosely identified same-flavor opposite-sign lepton with 76 < Mu < 106 GeV/c 2 . 
To maintain the projective geometry of the calorimeter, all events are required to have 
\z vtx \ < 60 cm. 

An inclusive £77 subsample is made by requiring each event to further contain a high pr, 
isolated, hadronically decaying r lepton candidate with pJp > 15 GeV/c [11] and \rj\ < 1.0. 
A 17 candidate is identified as a calorimeter cluster which satisfies the following require- 
ments [12]: (i) not identified as an e or /i; (ii) one or three tracks with pr > 1 GeV/c 
in a 10° cone around the calorimeter cluster center; (iii) the scalar sum of the p T of all 
tracks in AR = 0.4 around the cluster center, excluding those in the 10° cone, less than 
1 GeV/c; (iv) fewer than three ir° — >■ 77 candidates identified in the CES; (v) more than 
4 GeV of E T measured in the calorimeter; (vi) 0.5 < E T /pTp < 2.0 (1.5) for one track 
(three tracks); (vii) the cluster width of the calorimeter cluster in rj-cj) space less than 

8 



0.11 (0.13) - 0.025 (0.034) x E T [GeV]/100 for one track (three tracks); and (viii) the 
invariant mass reconstructed from tracks and 7r°'s less than 1.8 GeV/c 2 . The charge of the 
r h object is defined as the sum of the track charges, and is required to have unit magnitude 
and have the opposite sign (OS) of the I candidate. A total of 642 events pass the above 
requirements; 16 of these have two or more jets (identified using a fixed cone algorithm 
with AR = 0.4 [13]) with E T > 15 GeV and \rj\ < 2.4. Note that, as expected, the four 
lr h + jets candidate events which were found in the search for ti — > (W + b)(W'b) [12] pass 
the kinematic requirements for this search. 

The dominant backgrounds come from Z j^*[-^> T + r~)+jets, ti, diboson (W+W - , W ± Z, 
ZZ) production, and fake £r h combinations from VT+jets and QCD events. Monte Carlo 
(MC) programs with CTEQ4L parton distribution functions (PDFs) [14] and a detector 
simulation are used to estimate the background rates by simulating the kinematics of 
Z/Y, W, ti, and diboson events. All SM processes except W/Z+jets events are gener- 
ated using isajet [15]; VECBOS [16] is used for vector boson plus jets production and decay, 
followed by herwig [17] for the fragmentation and hadronization of the quarks and gluons. 
The cross sections for Z/Y, ti and W + W~ production are normalized to the CDF measure- 
ments [18-21] and next-to-leading order (NLO) calculations for W ± Z and ZZ production 
are used [22,23]. The number of fake events from QCD is estimated from the data and 
assumes that the number of OS events, after subtracting off the non-fake contribution, is 
identical to the number of like-sign (LS) events observed in the data as expected from QCD 
sources i.e., N<$ D = N^ ta - N^ c . 

The final optimized data selection requirements are based on simulated t\ti production, 
using isajet [15] and the CDF detector simulation, background expectations, and a control 
sample. See Fig. 2. To reduce the number of W+jets events we require Mt(£,$t) < 
35 GeV/c 2 where M T (£,fk) is the transverse mass of the I and the event fhr, defined 
as M T (£, fir) = ^2 p^fir(l — cos0£^), and where 0^ is the azimuthal angle difference 
between the t and the fir- To reduce the QCD backgrounds we require Sp T (£, Th,fir) = 
Pt + Pt + fir > 75 GeV/c. A control sample of £r h +0 jet events with similar kinematic 

9 



g 

> 



60 



50 



40 



30 



20 



10 











1 



CDF Run 1, 106 pb" 1 
• Data (OS Qj h +0 jets) 

□ z-wV 

□ QCD 



6 7 
N track (x) 



FIG. 1. The number of charged tracks in each candidate for the opposite-sign (OS) £rh+0 
jet control sample. The data are compared to the MC expectation (all background histograms are 
summed) which is dominated by real r^'s from Z — > t + t~ production. 



requirements (M T (t,fir) < 25 GeV/c 2 , + $r\ > 25 GeV/c) is selected to show that 
the backgrounds are well modeled, dominated by real Z — ► t + t~ production, and for later 
use in the acceptance calculations. Figure 1 shows the charged track multiplicity of the 
hadronic tau decays (removing the 1 and 3-prong requirements) for this sample and shows 
good agreement with background expectations. 

A comparison of the OS £rh+ > 2 jet data and background estimation is shown in Fig. 2 
before the final Mt(£,^jt) and £px(A T h, Pr) cuts. A breakdown of the backgrounds and 
data is given in Table I. The backgrounds appear well modeled. A total of 3. 21^3 events are 
predicted from all SM sources, dominated by Z{— > t + t~ )+jets production. No candidate 
events pass the final titi selection criteria, which is unusual but expected in roughly 3% of 
experiments when taking into account the statistical and systematic uncertainties. 

In order to set limits on t\ti production and decay, the acceptances and efficiencies are 
normalized to the rate of Z{— > t + t~)+0 jet decays using the following relation: 
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FIG. 2. The final data selection criteria for the OS £rh+ > 2 jet sample. The arrows show 
the final event selection requirements. The assumed stop mass is 100 GeV/c 2 . The quantities 
t/j, ftr) and Mt(1,$t) are defined in the text. 
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FIG. 3. The 95% C.L. upper limit on cross section for t\ti production compared to the NLO 
calculations 



'fail - r+r-bb) = ( ™^-J™^ \ • R Acc ■ R Trtg ■ a z ■ Br(Z -> t+t~) (1) 

where ^^Qp and (N^ bs and N^P) are the number of candidate events observed in the 

data and expected background in the > 2 jet/iiti (0 jet/Z) selections, Ra cc is the ratio of the 
Z to t\t\ acceptances and Rthq is the ratio of the trigger efficiencies. The primary advantage 
of this approach is that potential systematic uncertainties in the estimate of identification 
and isolation efficiencies are reduced in the ratio of t-\t\ to Z production. 

The 95% confidence level (C.L.) limits on a(titi — > T + r~bb) in the e, /i and combined 
channels are found using Eq. (1) and come from a Bayesian integration of the likelihood as 
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a function of the cross section, integrating over the correlated and uncorrected systematic 
uncertainties on the expected signal with a flat prior. The Racc term is a function of 
the ti mass and varies in the range 0.34 < R e Acc < 2.15 (0.35 < R\ cc < 1.87) for the 
e (/i) channel over the range 70 < < 130 GeV/c 2 . The Rrrig term varies between 
0.95 < R-Trig < 0-97 (0.99 < R% rig < 1.00) for the e (/i) channel with an uncertainty of 
about 1%. (The acceptance and trigger efficiencies for the Z control sample for this analysis 
are 1.19% (0.69%) and 74.5% (83.0%) for the e (/i) channel respectively.) Assuming lepton 
universality gives a z ■ Bi(Z -> r + r~) = a z ■ Bt(Z ->• = 231 ± 12 (stat+sys) pb [24]. 

The dominant uncertainty is due to the statistical uncertainty in N^ bs — Ajp and is 17.0% 
(24.9%) [25]. Additional uncertainty comes from our estimation of Ra cc which is dominated 
by the variation in the titi acceptance from choices of the QCD renormalization scale Q 2 , 
PDFs, amount of gluon radiation, the jet energy scale and the statistical uncertainty in the 
MC samples [26]. The total uncorrelated uncertainties vary between 17.1 and 17.7% (25.1% 
and 25.4%), and the total correlated uncertainties vary between 9.3 and 14.1%. 

Figure 3 shows the final 95% C.L. upper limits on the cross section times Br for the e, \i 
and combined channels, along with the NLO prediction of the production cross sections [27]. 
The 95% C.L. lower limits on Mr are 110 and 75 GeV/c 2 for the e and fi channels, respec- 
tively, where we have assumed Br = 1 for simplicity. Combining the two results yields a 
limit of 122 GeV/c 2 . Since our analysis does not distinguish the quark flavors in jet recon- 
struction, these results are equally valid for any X' 33k coupling. These results substantially 
improve on the currently most stringent mass limit which comes from the ALEPH exper- 
iment [28] which excludes t x masses below 93 GeV/c 2 using e + e~ — > t]ii — > t + t~+ 2 jets 
topology with an assumption of \' 33k ^ (k — 1, 2 or 3). 

In conclusion, we have searched for titi production using 106 pb -1 data in pp collisions 
at y/s = 1.8 TeV. We have examined the £r h + > 2 jet final state within an ^ SUSY scenario 
in which each t\ decays to a r lepton and a b quark via non-zero A333 or £3 couplings. No 
titi event candidates pass our selection criteria and we have set a 95% C.L. lower limit on 
the h mass at 122 GeV/c 2 for Br = 1. 
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Sample 


tt 


Diboson 


W + jets 


Z/ 7 * -»• r+r- 


QCD 


Tot 


Nobs 


OS It 


1.2±0.3 


2.3±0.8 


101±6 


225±9 


301±18 


631±21 


642 


tT h + > 2 jets 


1.0±0.2 


0.4±0.1 


3.4±0.4 


7.7±0.5 


8±3 


21±3 


16 


M T {l,$r) < 35 GeV/c 2 


0.15±0.07 


0.14±0.06 


0.5±0.2 


6.0±0.4 


8±3 


15±3 


10 


Z PT (£,T h ,fir) > 75 GeV/c 


0.15±0.07 


0.08±0.03 


0.24=0.1 


2.8±0.3 


n+1.4 
u -0 


O 9 + 1.4 

°- z -0.3 






TABLE I. Summary of the number of OS events in the data and expectations for the background 



sources as each selection requirement is applied. 
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